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We investigate energetically favorable configurations of point disclinations in nematic 
films having a bump geometry. Gradual expansion in the bump width A gives rise to a 
sudden shift in the stable position of the disclinations from the top to the skirt of the 
bump. The positional shift observed across a threshold A th obeys a power law function 
of | A — A th | , indicating a new class of continuous phase transition that governs the 
defect configuration in curved nematic films. 
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1. Introduction 

Optimal control of molecular alignment is crucial for applications of liquid crystal 
films. Thus far, many sophisticated techniques have been developed for obt aining; 
homogeneous alignment in a d esired direction: mechanical rubbing trea tment J ' \ 
surface-polarization effect surface coating ^ and photoalignment I ^M-M I are 
only a few examples to mention. When employing these techniques, topological de- 
fects are regarded as unwanted matters and thus they are eliminated artificially. 
Under certain conditions, however, topological defects can have beneficial effe cts on 
the alignment control of liquid crystal molecules 1 1 1 1 1 2 | 1 3 | 14 | 1 5 | 1 6 | 1 1 8 19 20 as well 
as of guest particles embedded in a liquid crystal medium | 21 | 22 | 23 | 24 | 25 | 26 | 2 7|28| 
The positive effects result largely from topological constraints for the molecular 
orientational field (= director field) that is confined in a closed curved substrate. 
In fact, the director field over closed curved surfaces strongly depend on the posi- 
tions and charges of topological defects distributed on the surfaces, which implies 
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the controllability of the director field through the defect manipulation. Such the 
utilities of topological defects may be enhanced when combined with the geometric 
potential effect. It is known that topological defects in curved liquid crystal films 
experience a surface-curvature-indu ced potential, attr active or repulsive depending 
on the local geometry of the surface I 29 | 30 | 31 | 32 | 33 | 34 ]ja| Therefore, a better under- 
standing of the correlation between curvature variation of the underlying surface 
and favourable configurations of topological defects will shed new light on molecular 
alignment control of liquid crystal films. 

A Gaussian bump is one of the exemplary curved surfaces for studying the cor- 
relation mentioned above. It has a continuous rotation symmetry around a vertical 
axis, showing a positive (negative) Gaussian curvature close to (far from) the top 
of the bump [see Fig. [T^a)]. The coexistence of a positively curved region and a 
negatively one on the same surface allows a simultaneous observation of the curva- 
ture effects at the two distinct regions with oppositely-signed Gaussian curvature. 
Earlier work suggested that when a nematic liquid crystal film is placed on a Gaus- 
sian bump, a point disclinatio n with +1 charge is trapped to the top of the bump 
due to the curvature potential I 32 | 33 [ Contrariwise, we have found in the previous 
work that point disclina tion s are attracted to the inflection point of the bump re- 
gardless of their charges S^l These two results appear controversial; a possible clue 
to resolve the controversy is the difference in the assumed bump widths between 
the two existing studies. In the former study, the bump width is sufficiently small, 
so that the bump is covered almost entirely by a single disclination as de picted in 
Fig. [He). This situation is expressed by saying that the coherence length EsMQJ £ ( 
of the nematic order fluctuation exceeds the scale of the bump width (see Sj2] for 
the definition of £). In contrast, the latter study was based on the assumption that 
the bump width is larger enough than the coherence length and the system can 
involve many small disclinations [see Fig.fljb)]. It is thus reasonable to conjecture 
that the ratio of the bump width to the coherence length plays a decisive role in 
determining the stable position of disclinations on a Gaussian bump. Furthermore, 
the two contrasting results imply the existence of a threshold ratio across which 
the stable position of disclinations alters. 

In the present Letter, we prove our conjecture by showing a sudden shift in 
the stable position of a +1 disclination with increasing the bump width. Once the 
bump width exceeds a threshold value, the stable position is suddenly switched 
from the top to the inflection point of the bump. The positional shift is described 
by a power law function of the bump width, as is analogous to the second-order 
phase transitions of various physical systems whose properties are characterized 
by power-law behaviors. The present result thus indicates a novel kind of phase 
transition that is peculiar to curved nematic films. 

a Similar cur vatu re effects were discussed in the literature of quantum transport | 35 | 36 | an( j c ] ass i ca i 
spin models on curved surfaces. 
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(a) (b) . a (c) 




Fig. 1. (a) Sketch of a Gaussian bump characterized by the width A and the hight v. The region 
inside (outside) the dashed circle with radius A possesses a positive (negative) Gaussian curvature, 
(b, c) Disclination embedded in the Gaussian bump. The coherence length £ of the nematic order 
fluctuation measures the spatial extent of dcsclinations (hatched regions) in which the director 
field n(R) is arranged in a radial manner. 



2. Method 

We consider point disclinations in a nematic liquid crystal membrane having a 
Gaussian bump geometry. The bump is described by a position vector R(r, <f)) — 
rcos 4>i+r sin0 i 7+wexp(— r 2 /2A 2 )fc, where k are right-handed orthogonal basis 
vectors [see Fig.QJa)]. The constant A, which we call the bump width, determines 
the inflection point of the bump at which (d 2 R/dr 2 ) ■ k — 0. A dimensionless 
parameter a = v/A quantifies the convexity of the bump; a = for the flat plane, 
and a 1 for a spiky bump. We fix a = 0.1 throughout the paper, as a result of 
which spatial variation of the surface curvature is smooth and sufficiently slight to 
be accessible in experiments. 

Our aim is to evaluate the stable position of a disclination constrained on the 
Gaussian bump. It is solved by calculating the elastic free energy F<j of the director 
field n(_R) around the disclination core, since takes a minimum at the stable 
position to be evaluated. We assign the core of a +1 disclination at R c = R(r c , <f> c ) 
with arbitrary (j> c , and assume that the disclination is endowed with the coherence 
length £ of the nematic order fluctuation. The latter assumption ensures a radial 
pattern of n(_R) within a circular region A having the radius £ centered at the core 
[i.e., the hatched regions depicted in Figs. QIb) andQJc)] El- The radial pattern of 
n(R) is gradually distorted with increasing distance from the core; in fact, two-point 
correlation in the director orientation disappears if the points are well separated by 
a distance larger than £. A typical value of £ is on the order of tens microme ters as 
can be estimated from polarization microscope images of liquid crystal films 1 41 1 42 | 

The contin uum elastic theory states that associated with the disclination 
rea ds as™] 

^ d = hJdA^{D ■ nf + ^(D x n) 2 + ■ [(n • D)n - n(D ■ n)] J(I) 



b Precisely speaking, the area A slightly deviates from an exact circle due to the finite curvature 
of the surface. In fact, we define the perimeter of A by the locus of points whose distances from 
the core are equal to £ in the sense of geodesic distance. Therefore, the perimeter is distorted a 
bit from the exact circle. 
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Fig. 2. (a) Disclination energy in deformed MBBA films as a function of the core position r c . 
(b) Sudden shift in the stable position r™ ln as a function of the bump width A. The threshold 
A th is estimated to be A th /£ = 0.95. Inset: The same data in the log-log plot, showing a power 
law r™ in tx | A - A th | 8 with the exponent /S ~ 0.2. 



where D is a vector operator whose components are the covariant derivative 
n the curved surface. The constant h is the thickness of the liquid crystal film; 
Kx, K 3 and K 2 4 are the elastic constants associated with splay, bend and saddle- 
splay distortions of the field n, respectively The integration range in Eq. ([TJ is 
restricted to the finite circular region A, despite of slow decay in the elastic energy 
density with distance from the core; we will revisit this issue in $4] To avoid the 
divergence in at the core, we introduce an inner cut-off A corc around the core 
according to the formula A COIC = irKi/2kB(T c — T), in which T c is the nematic- 
isotropic transition temperature and T is the temperature of the system S3. 

In actual calculations, we set £ = 20/im and h = 4 /im by referring to 
the experimental observation BOHM ^he values of K\ and K3 are those of 4- 
methoxybenzylidene-4 / -butylaniline (MBBA): K\ — 6 pN, K3 — 7.5 pN at room 
temperature Since the corresponding value of K24 was not yet measured, we 
set K24 = 4.9 pN on the basis of the relation K24 = {K\ + K 2 )/2 that involves the 
twist coefficient K 2 = 3.8 pN for MBBA S3. 



3. Results 

Figure [2ja) shows the disclination energy as a function of the core position r c . 
The normalized bump width, A/£, ranges from 0.8 to 2.0. For A/£ = 2.0, the curve 
of Fd is convex downward, showing a minimum at r c ~ A. The energy scale of 
the potential depth corresponds to U-qT at a temperature of T = 10 4 K, which is 
by far larger than the thermal excitation energy at room temperature. Hence, this 
downward peak implies that a disclination is trapped to the annulus region having 
the radius of r = A surrounding the bump top [see Fig.^a)]. We call this annular 
distribution of disclinations by the "halo" phase, according to the discussion in Ref. 
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What deserves attention is the A-driven transition of the downward peak profile. 
It follows from Fig. [2ja) that both the position and magnitude of the peak are 
insensitive to the decrease in A/£ until it reaches a specific value of A/£ = 0.95. 
A drastic change takes place at A/£ = 0.95; the downward peak in Fi is replaced 
by a plateau that extends from r c = to r c ~ 0.5A. This plateau means that the 
disclination can move freely within a finite circular region centered at the bump 
top. It is also found that a further decrease in A/£ gives rise to a new downward 
peak in Fi at r c = 0. That is, when the bump width A is comparable to (or smaller 
than) the coherence length £, then the disclination is trapped at the bump t op [see 
Fig. EJb)]; this result is in complete agreement with the earlier observation I 32 | 33 ] 
We call this defect configuration by the "convergence" phase. 

We have seen that the peak of Fi appearing at r c = r"" n is gradually shifted 
from the inflection point (r™ 111 ~ A) to the top (r™ m = 0) of the bump with 
increasing A. Figure 2(b) shows the A-dependence of the normalized peak position 
r" lln /A. For a large bump width (A/£ ;§> 1), Fa takes the minimum at r" lm /A ~ 1 
implying the halo phase. The disclination halo shrinks with decreasing A, and 
finally it converges at the bump top when A = A th = 0.95£. These results indicate 
the continuous phase transition (i.e., of the second order) of the stable disclination 
configuration on the Gaussian bump; the system undergoes the transition from the 
halo phase to the convergence phase across the threshold A th /£. This is the main 
finding of the present Letter. 

It is interesting to note that the A-dependence of r™ ln in the halo phase is 
described by a power law of the form 

r™ n oc | A- A th f . (2) 

The inset of Fig. 2(b) evidences the power-law behavior of Eq. ([2]) with the exponent 
j3 ~ 0.2, which holds over the two orders of the lateral axis. We thus believe that 
the present result is a precursor of an untouched configurational phase transition 
peculiar to curved liquid crystal systems, though further studies are needed to 
establish its universality and physical interpretation. 

4. Discussions 

The transition of the stable defect configuration we have demonstrated is a conse- 
quence of the following two Fi variations caused by reducing A: a slight Fi-increase 
to the right of the arrow P marked in Fig. 2(a), and a rapid Fi-decrease to the left 
of P. In this section, we give an outline of the physical mechanisms of the two Fi 
variations, while a detailed formulation that describes them will be given elsewhere. 

The rapid decrease in Fi to the left of P originates from a suppression of the 
splay deformation of the director field n(R) around the bump top. Suppose that 
a +1 disclination is placed near the top of a broad bump as depicted in Fig. [TJc). 
Since £ and a — v/A are constants, a reduction in A leads to a lowering of the bump 
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(a) Halo phase (A>£) (b) Convergence phase (A<Q 




Fig. 3. Distribution patterns of disclinations in the halo phase (a) and convergence phase (b). In 
either phase, disclinations are trapped in the dark region. 



hight v remaining the disclination area A be fixed. As a result, the disclination is 
allowed to cover the whole bump. In the latter situation, adjacent directors (thick 
line segments in the plots) distributed close to the outer edge of A tend to be 
parallel to each other. Eventually, the first integrated term in Eq. fTJ), oc (D ■ n) 2 , 
which represents the contribution from the splay deformation of n, is suppressed 
by reducing A. This is the reason why at r c = drops off when A decreases. 

A parallel argument to the above accounts for the slight increase in to the 
right of P in Fig. [2ja). Suppose a +1 disclination on the skirt (i.e., at r c ~ A) of 
the bump. If the disclination area A intersects the circle r = A around the bump 
top, then it covers a positively-curved region (r < A) and a negatively-curved one 
(r > A). At the negatively curved region, a reduction in A yields an enhancement 
of the splay deformation which leads to the slight increase in F^ at r ~ A. A 
mathematical proof of the splay enhancement has been given in Ref. 1221. 

The data shown in Figs. [2](a) and [2](b) are derived from a specific Gaussian 
bump shape characterized by the parameter a = v/A = 0.1 (see Fig. Ufa)). We also 
confirmed that, when 0.05 < a < 0.5, the system exhibits the same configurational 
transition across A th . A slight increase in A th with a was observed within the 
range of a, reflecting further splay suppressions and enhancements in disclinations 
whose cores locate at r c = and r c = A, respectively. In contrast, the value of f3 is 
insensitive to the change in a. 

It should be stated that our discussion is based on the model given by Eq. (JI]), in 
which contributions from the outer region than the circle A are ignored. In actual 
nematic layers, a +1 disclination has a profile of the elastic energy density that 
decays logarithmically with distance from the core. Hence, contributions far from 
the core may be relevant to the stable position of the disclination, thus should be 
taken in account in order to obtain a quantit ativ ely accurate value of A th . The 
geometric potential approach developed in Ref. ^Hl will be a clue to solve this issue. 

Before closing the Letter, we suggest the following system could allow experi- 
mental tests of our theoretical predictions; it is a thin liquid crystal film coating 
a randomly corrugated substrate. The film involves many disclinations and bumps 
with various sizes; therefore, we will observe simultaneously the two effects, halo 
and convergence patterns depicted in Fig. [3J by tuning the surface corrugation. We 
also mention that such the corrugated substrate can be synthesized by polymer 
adhesion with patterned surfaces as demonstrated in Ref. ^4 
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5. Conclusion 

We have found a continuous transition of defect configurations in nematic films 
confined on a Gaussian bump. Depending on the ratio A/£, the system exhibits 
either of two defect distribution patterns: the halo phase (A/£ > 1) and convergence 
phase (A/£ < 1). The stable position r™ m of defects turned out to obey the power 
law r™ m oc | A — A th |^ with the exponent (3 ~ 0.2 under the present condition. Our 
findings suggest the possibility of defect manipulation via artificial deformation of 
liquid crystal films. 
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